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SUMMARY 

Owing to the quantity of material to be presented, this 
report is divided into several volumes. Volume 1 introduced the 
series and presented and discussed the results for Shot 10. 
Volume 2 presents and discusses the results for Shot 9. 
Subsequent volumes will present and discuss the results for 
Shots 8 and 1l. The method of analysis is common to all four 
experiments and is described in detail in Volume 1 only. 

So that the results from the four experiments may be 
easily compared, they have been scaled to remove the effects 
of varying atmospheric conditions. (Results are scaled toa 
1 kg charge weight and a standard atmosphere of dry air at 
15°C at sea level.) For the most part, only scaled results 
are presented. Exceptions include some derived pressure-time 
histories, which may be compared to actual gauge measurements 
made in the experiment. 

Results are presented in SI units, even though the experi- 
ments were originally laid out in British units. Only distance 
and time measurements are affected, however, as velocity 
density, and pressure results are presented as dimensionless 
ratios. A distance units conversion scale is included on page 3 
to convert between SI units (meters scaled to a 1 kg charge) 
and British units (feet scaled to a 1 1b charge), plus a time 
scale factor and scale factors to convert pressure ratios 


to both British and SI pressure units. Scale factors 


ees 


which may be used to compute the distance and time values 
actually observed under the ambient conditions of each shot 


are also provided. Dimensional pressure units are used for 


the results presented at gauge locations. 
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CHAPTER 1, SHOT 9 ANALYSIS 


1.1 Introduction 

This is the second volume in a series which presents the 
particle trajectory analysis results from four experiments 
(Dipole West Shots 8, 9, 10 and 11) carried out to obtain 
information on the interaction of spherical blast waves with 
real and idealreflecting surfaces. A general description of 
the project can be found in Volume 1. The results presented 
in this volume are for Shot 9, in which the same charge 
configuration was used as in Shot 10 (reported in volume l), 
but which was carried out over a smoother ground surface. In 
each experiment, photogrammetrical studies were made of the 
shock frents (refractive image analysis, RIA), and of the 
motions of smoke puff particle tracers (particle trajectory 
analysis, PTA). The refractive image analysis results were 
reported by Dewey et al. (1975) and results of the particle 
trajectory analysis are presented in this report. The method 
of particle trajectory analysis, common to all four shots is 


described in detail in Volume 1 only. 


1.2. Description ‘of Shot 9 
Dipole West Shot 9 was fired on October 22nd, 1973 by 


the Ballistics Research Laboratories at the Defence Research 


Establishment Suffield, in Alberta, Canada. Two 1080 lb 


(491 kg) spheres of Pentolite were detonated simultaneously, 
to within 5 microseconds, at nominal charge heights of 15 and 


45ft (4.6 and 13.7m) over a relatively smooth ground surface. 


: 
= | 
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Particle trajectory data were gathered by photographing 
the movement of smoke puffs formed in a vertical plane running 
out from ground zero at 6.7° south of west. A WF5 camera 
Operating at about 3400 frames per second was positioned 
30ft (9.2m) above ground level at a position 600ft (183m) 
due south of ground zero (GZ), the point on the ground vertically 
beneath the charges. 

Table 1 gives the field survey data for the event, and 
Figure 1 shows a plan view of the layout. The dashed line 
represents the approximate line of sight of the WF5 camera. 
Figure 2 shows the field of view of this camera. 

The smoke puff grid was made up of 9 columns of 12 puffs 
each, hung vertically on strings. The vertical spacing of 
puffs was 5ft, beginning 3ft above ground level and ending 
at a height of 58ft. The horizontal spacing of the columns 
of puffs was 10, 7 or 5ft, depending on the distance from 
ground zero, beginning at about 25ft and ending at about 85ft 
from GZ. Of the possible 108 smoke puffs, 106 detonated 
successfully. A good film record was obtained, except that 
the detonation zero timing mark was not recorded, which means 
that the results cannot be related to the detonation pulse 


to the charges with an accuracy better than +0.15 ms. 
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This report describes the analysis of the smoke puff 
data collected for Shot 9, and presents and discusses some 


of the results of that analysis. 


1.3. Camera calibration and data reduction 

The calculated camera position coordinates and orien- 
tation angles for Shot 9 are presented in Table 2, together 
with the positions of photomarkers transformed from one frame 
of the film just before detonation to an object plane defined 
as passing through ground zero and being normal to the camera 
orientation axis. The differences ("Shifts") between the 
object plane positions of the transformed calibration points 
and their positions computed from the field survey data are 
given in Table 2. The object plane positions of the calibra- 
tion points computed in these two ways are also shown in 
Figure 3. 

The camera calibration procedure, described in detail 
in Volume 1, ensured that selected photomarker images (Pl to 
P5) transformed to the object plane in a way which matched 
them exactly to the positions computed using the survey data. 
These reference photomarkers for Shot 9 are indicated in 
Figure 3 using large circles: namely, Pl = Wl, P2 = W3, 
and P3 = 300W2. The separation distance between P4 = P3 = 
300 WZ and PS = W2 was also used as a calibration parameter. 


The probable reasons for the shifts seen for photomarkers VP1A 


1] 


and VP1B was discussed in Volume 1. 

The image positions of two reference photomarkers 
(VP3B and 300W2) and all smoke puffs were measured frame-by- 
frame over a time interval corresponding to the approximate 
duration of the positive phase of the blast waves (film frames 
9 to 200), and were transformed to distances in the object 
Plane by matching the reference marker positions to their 
positions transformed from the calibration frame. These data 
were again transformed from the object plane to the smoke 
puff plane which was assumed to pass through "corrected" 
ground zero; to be vertical, and to run 6.7° south of west 
from GZ. 

The x-y coordinate system in the smoke puff plane 
was the same for Shot 9 as for Shot 10 except that the corrected 
value for ground zero was displaced 1.0 ft from the surveyed 
ground zero, in a direction approximately 46° south of east. 
The corrected ground zero was defined to have the same 
elevation as the surveyed ground zero, but was located directly 
under the midway point between the two charge centers. As 
for Shot 10, all data in the output plane are plotted with 
the x coordinate reflected, i.e. with positive values of x 
to the right hand side, as if the smoke grid had run to the 
right of the charges rather than to the left as seen in the 
film images. 

A time was assigned to each film frame using the 1 ms 
timing marks placed on the film during its exposure. The 


film timing method was described in Volume 1, and the complete 


set of film timing data used for Shot 9 is provided in 


Table 3. In the absence of a zero time pulse on the film, 
zero time was determined using the static zero calibration 
distance for the camera plus one-half a frame length. Because 
the static zero distance is at most a frame length less than 
the actual zero distance measured from any moving film, the 
error in frame times can be at most one-half frame (about 

0.15 ms). 

Figure 4 shows the positions of the 106 detonated smoke 
puffs at a time prior to the detonation of the two charges. 
These positions are in the plane of the charges and the 
smoke puff grid, as described above. The smoke puff plane 
was not exactly parallel to the camera image and object 
planes (Figures 2 and 3), and various geometrical corrections 
were applied to make the transformation between them. The 
puffs enclosed in parentheses were not visible in the earlier 
film frames, but were seen later when they were illuminated 
by the light of the fireball. Charge positions in the figures 
are plotted as if they were positioned exactly above the 
corrected ground zero origin. The data shown in Figure 4 


have not been scaled. 


1.4 Data scaling and trajectory fitting 


The position-time histories of individual smoke puffs were 
extracted from the frame-by-frame positions of the smoke 


puff grid, and then scaled to standard atmospheric conditions 
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and charge weight. 


point in the analysis. 


edited, and then smoothed by fitting polynomial functions. 
Particle trajectory data were scaled by dividing all 
= Y(wP)AW P) 
Ss (WPOIAWSP) and 


C/(C)S) , where C is 


distances by Sachs scaling factor 


multiplying all times by the factor 


A change to SI units was made at this 


The resulting trajectories were 


the ambient sound speed computed for Shot 9. 


compute C and S, and define the scaled event, 


below with the computed values of C and S. 


Ambient temperature, 
Ambient pressure, 
Relative humidity, 
Computed vapour pressure, 


Ambient sound speed, 


Charge weight, 


Sachs scaling factor, 


Standard charge weight, 
Standard pressure, 
Standard temperature 


Standard sound speed, 
(dry air) 


PT = £4.17 °C 


P = 93.02 kPa 


RH = 55.0% 


VP = 0.89 kpa 


C = 340.469 m/s 


W = 489.9 kg 


6S = 8.1111 

Wo = 1.0 kg 

Ps = 101.325 kpa 
— ° 

7. = 15 °¢ 

Cy = 340.292 m/s 
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Data used to 


are listed 


(57.5 “F) 


(13.491 PSL) 


(6.7 mm Hg) 


(1117 £t/s) 


(1080 lbs) 


(2.2 UBS) 
(14.7 PST) 
(59 °F) 


(1116 £t/s) 


The results presented in this report therefore apply 
to a scaled event which is the detonation of two 1 kg charges 
in a standard atmosphere. The scaled heights of burst for 
Shot 9 were 0.571 m and 1.707 m, and the charge separation 
divided by two, scaled, was 0.568 m. These figures may be 
compared to the scaled charge height and half-separation 
distances for Shot 10 which were 0.563 and 0.575 m respectively. 

Figure 5 shows the scaled particle trajectory data for 
Shot 9 in the smoke puff plane with positions measured horizontally 
ane vertically from corrected ground zero. Approximately 
9350 puff positions are represented. As represented, the raw 
trajectory data havenot been smoothed. 

The raw particle trajectory data were edited to remove 


obvious data processing errors, such as a single point 


widely displaced from its trajectory for one or two frames. 
The trajectory of each puff in turn was then smoothed by 

least squares fitting simple polynomial expressions separately 
to both the x and y coordinate data, these being discrete 
functions of frame time. The adequacy of each fit was 
determined by examining on the same graphical output, plots 

of both the raw trajectory data and the fitted curve. For 
Shot 9 this meant examining and adjusting 212 such plots, at 


least two or three times each. 


For a given puff, the first step in fitting the raw 


trajectory data was to set the time of arrival of the shock 
front first hitting the puff. The data at subsequent times 
were fitted with polynomial functions, as described in 
Volume 1, paragraph 2.5. The first derivatives of the fitted 
functions were also calculated at a series of times for use 


in later calculations of particle velocity. 


1.5  Regionalization and shock strength calculations 


Five regions were defined in the smoke puff plane on 
the basis of the shock front which first struck the puffs in 
a particular region. These are shown in Figure 6. The 
regions were bounded by the triple point trajectories measured 
using refractive image analysis (Dewey et al., 1975). Regions 
1 and 2 are those in which the smoke puffs were first hit by 
a spherical primary shock front, and regions 3, 4, and 5 
are those in which the puffs were first hit by a Mach stem. 

In each of the five regions, the shock trajectory data 
obtained from the first movement of the smoke puffs were fitted 


to a function of the form 
EUG) =A Bere log (i + ec) | 


where r is the shock radius, t is the time after detonation, 


and A, B, and C are the fitted coefficients. The shock 


velocities were calculated by differentiating this function. 


The peak particle velocity, V 


3 i peak density, De , and peak 
hydrostatic overpressure, i , aS functions of shock radius 
in each of the five regions, were calculated from the shock 
velocity using extensions of the Rankine-Hugoniot equation. 


Details of the shock radius calculations etc. are described 


in Volume 1, paragraph 2.6. 


1.6 Particle velocity calculations 


Particle velocities were computed using the methods 


described in Volume 1, paragraph 2.7. 


1.7 Density and hydrostatic overpressure calculations 


Densities and hydrostatic overpressures in the smoke puff 
plane were calculated by the method described in Volume l, 
paragraph 2.8. Results in both cases represent average 


values over cells defined by four adjacent smoke puffs. 


1.8 Surface representation 


Surfaces were fitted to the times of shock front arrival 


and to the fields of particle velocity, density and hydrostatic 


overpressure at a sequence of times. All data were interpolated 
onto a common regular Euleurian grid. Fields of dynamic 
pressure were computed from surface-interpolated particle 
velocity and density results. Contour plots were generated for 
all surfaces at selected times, and time histories computed at 
several fixed locations. The methods used were identical to 


those described for Shot 10 in Volume 1, Chapter 3. 
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CHAPTER 2. SHOT 9 RESULTS 


2.1 Times of shock front arrival 

The measured initial puff positions, the times of first 
shock arrival, and the peak particle velocities obtained by 
differentiating the functions fitted to the particle 
trajectories are presented in Table 4. Puff position is 
Given relative to corrected ground zero as origin with 
horizontal and vertical axes. Puff position and the time 
of arrival of the first shock are given both as observed and 
scaled. Particle velocities listed are derivatives of the 
fitted puff trajectories at the times of shock arrival, and 
are expressed in Mach units. Expressed this way, the particle 
velocities are the same scaled as unscaled. Also listed are 
the initial radial puff positions (scaled) and region codes. 

Shock front data determined from the first movement of 
the smoke puffs, i.e. calculated from the time-of-arrival 
data in Table 4, are listed in Tables 5.1 - 5.5. Each 
table corresponds to one of the 5 regions used. Listed are 
the observed and fitted unscaled shock trajectory data, the 
scaled fitted shock trajectory data, and the computed shock 
velocities and peak parameters associated with shock 


strength: peak hydrostatic overpressure in atmospheres and 


in kilopascals, peak particle velocities in Mach units, and 


peak density ratios. Given as ratios, these peak parameters 
are the same scaled as unscaled. Pressure given in kilo- 
pascals in the tables refers to the unscaled (observed) case 
only. 

The shock front radius versus time data derived using 
particle trajectory analysis (PTA) are also shown in Figures 
7.1 - 7.3 for the two primary fronts, the two Mach stems at 
the interaction plane, and the ground Mach stem, respectively. 
They are compared to corresponding data derived from refractive 
image analysis (RIA) reported by Dewey et al (1975). The 
refractive image analysis results were obtained using photo- 
grammetry against a striped canvas backdrop and they describe 
the shock as it travelled in a direction almost diametrically 


opposite to the direction of the smoke puff grid. 


Z.2 SHOCK strengths 


Peak particle velocities calculated from shock front 
velocities are shown in Figure 8.1 - 8.3 for the primary 
fronts, interaction Mach stems, and the ground Mach stem. 

This method of determining peak particle velocities has 

been labelled method 1, and the data plotted correspond to 
those listed in Tables 5.1 - 5.5. The results in the figures 
are compared with those previously obtained using refractive 
image analysis (RIA). In the case of the primary shock fronts, 


results are also compared to those of Brode (1957) for TNT. 


In Volume 1 other methods of determining shock strengths 
in the various regions were described. It was demonstrated 
that method 1 was clearly the most accurate, and in the present 
volume shock strengths calculated using methods 2 and 3 are 
not reported. For this reason Figures 9, 10 and 11 and 


Table 6 do not appear in this volume. 


2.3 Particle velocity fields 

The calculated particle velocities in the plane of the 
smoke puffs are shown as vectors in Figures 12.1 through 12.7, 
for various times after the detonation. All times and 
positions are scaled to a 1 kg charge in a standard atmosphere. 
The particle velocity vectors represent the derivatives of 
the smoothed particle trajectories, and their magnitudes 
may be judged using the standard vector shown on each figure. 
All velocites are measured in Mach units, relative to the 
standard sound speed. Puffs not yet struck by a shock wave 
are represented by small circles (zero velocity). 

Numerical data corresponding to Figures 12.1 - 12.7 
are listed in Tables 7.1 through 7.6, along with scaled 
radial positions of the puffs, and region codes as defined 
in Figure 6. Conversion factors are given at the foot of each 


table, which may be used to convert the scaled data in the 


tables and figures back to their original unscaled values. 


2.4 Density and hydrostatic Overpressure fields 


Calculated average relative densities throughout the smoke 
puff plane are depicted graphically in Figures 13.1 - 13.4, for 
various times after the detonation. All time values are 
scaled. Cell positions are scaled and are given relative 
to the corrected ground zero as origin with horizontal and 
vertical axes. The calculated densities may be judged using 
the density shading scale shown on each figure. Density 
is given as a ratio, relative to ambient density. Cells 
not yet struck by a shock wave and cells in which the density 
has dropped to a value less than ambient density are shown 
blank. 

Corresponding numerical data are listed in Tables 8.1 - 
8.3 along with radial cell positions computed according to 
the regions defined previously. Numerical data describing the 
fields of hydrostatic overpressure are similarly listed in 
Tables 9.1 - 9.3. The pressure results for a given cell 
were obtained by multiplying the density results for that 
cell by a factor determined by the strength of the shock 
which first traversed the cell and which then remained constant, 


i.e. by assuming isentropic flow after the first shock. 


2:25 Times-of-arrival surface 
Figure 14 shows a perspective view of the surface fitted 
to the original smoke puff positions and the observed times 


of first shock front arrival, i.e., to the data listed in 


Table 4. The grid mesh size is 0.1 by 0.1 meters (scaled), 
about 2.5 feet square (unscaled), or about 4% that of the 
Original smoke puff grid. The charge positions are indicated 
on the vertical distance axis. 

The times-of-arrival surface is smooth enough to permit 
contouring, the contours in this case (isochrones) representing 
shock front shapes at different times, as shown in Figure 15, 
but the surface is not smooth enough to permit the calculation 
of gradient vectors which could be used to compute shock 
velocity vectors and shock strengths over the new grid. 

Two attempts were made to obtain contours of shock 
strength. In the first, the times-of-arrival surface was 
smoothed by least-squares fitting low-order, one-dimensional 
polynomial functions to the time~of-arrival data along each 
grid row and column separately, and computing the derivatives 
of the fitted functions to obtain the associated components 
of the surface gradient vectors. Shock velocity vectors were 
obtained from the time-of-arrival gradients, and from these 
peak particle velocities were computed. The peak particle 
velocity (shock strength) surface is shown in Figure 16. 

The contours of this surface (not shown) did not exhibit any 
discontinuities across the boundaries of the shock front 
regions, as they would if surfaces were fitted to the time of 


arrival in each region separately. 


bo 
bo 


The results of a second method used to compute shock 


strength contours are shown in Figure 17. These were obtained 
by interpolating shock radius at each value of peak particle 
velocity shown, for each shock front region in turn, using 

the peak particle velocity versus radius curves shown in 

Figure 8. Arcs of circles with these radii, centered on the 
appropriate points along the vertical charge axis, were then 
drawn in the regions to represent shock strength contours. 
These peak value contours are discontinuous across triple point 
locii and other region boundaries. As a result, some horizontal 
lines are crossed twice by the same contour or, in other words, 
identical shuck strengths can be found at two locations the 
same vertical distance from a reflecting surface, but at 


different radial distances from the vertical charge axis. 


2.6 Field surface contours 

Contours of equal particle velocity, density, hydro- 
static overpressure, and dynamic pressure in the blast 
waves were determined for a series of times, using surfaces 
fitted to the various measured data fields at those times. 
Sample results are shown in Figures18 through 21 at scaled 
times of 2.5 ms and 4.0 ms. The shock fronts shown in 
these figures are obtained from the time-of-arrival surface 


(as were these in Figure 15). Field contours such as those 


shown can be drawn for any scaled time between 0.5 ms and 7.0 ms. 


Sulieeareenerentwens 


It should be re-stated that all of these results were obtained 


from the photography of the smoke puffs only and do not rely 
on the results obtained using the refractive image analysis 


(Dewey et al., 1975). 


2.7 Time histories 

By mapping the physical properties of the blast waves 
at short time intervals it was possible to determine the time 
histories of these properties at any selected fixed position 
within the smoke puff grid. This was done at 12 fixed 
locations, three in the two primary regions and three in each 
of the three Mach stem regions, as shown in Figure 22. At 
each distance from the axis of the charges in the Mach stem 
regions, each of the time history stations is the same distance 
from either the interaction plane or the ground plane. Particle 
velocity time histories could be interpolated closest to the 
ground level because these were measured at puff locations, 
whereas the density and pressure data were measured at cell 
centers. 

Time histories of particle velocity, density, hydrostatic 


and dynamic overpressure at these locations are given in 


Figures23 to 26. 


Also plotted with the time histories are the interpolated 


values of the time of arrival of the first shock front at the 
stations. The height of this time-of-arrival line represents 
a peak value derived from the shock velocity analysis. 

Time histories for hydrostatic and dynamic pressure are 
also plotted in Figure 27.1 to 27.7 for stations at the 
nominal positions of field-mounted pressure gauges on the 
"60 foot gun barrel". The gauges on this gun barrel were 
mounted at nominal elevations of 10, 15, 20, 27, 30, 33, and 
40 feet. The time histories at these locations are given in 
unscaled units in order to facilitate comparisons with the 
gauge measurements. 

The dynamic pressures plotted in figures 26 and 27 are 
maximum values, computed using both the x and y components of 
particle velocity. Similar plots were made of the horizontal 
components Of dynamic pressure, but the differences were not 
significant since the y components of particle velocity at these 


locations were small. Other locations could have been chosen 


at which the y components would not have been insignificant. 


CHAPTER 3, DISCUSSION 


3.1 Particle trajectory analysis, Shot 9 


The methods used to analyze the smoke puff trajectories 
| on Shot 9 were identical to those used for Shot 10 and 
| described in detail in Volume 1 of this report. However, 
the results for Shot 10 clearly indicated the superiority 
of one of several methods of analyzing shock strength, and 
Only the results of the superior method were reported for 
Shot 9. 

The rate of smoke puff failure was greater for Shot 9 
than for Shot 10 (there were no failures for Shot 10), but 
the incompleteness of the Shot 9 particle trajectory grid 
did not greatly interfere with the analysis of the blast 
wave within the smoke puff region. The extra interpolation 
required for Shot 9, however, will undoubtedly correspond 
to some slight decrease in the accuracy of the calculated 
blast parameters compared with those for Shot 10. 

3.2 Primary shock strength of upper charge 

The refractive image analysis of the shock fronts 
described by Dewey et al, 1975 did not provide any information 
about the primary spherical shocks from the upper charges, and 
it was assumed that these charges had detonated satisfactorily. 
This assumption has now been validated for Shot 9 by the 


analysis of the particle trajectory time-of-arrival measurements. 


In Figure 7.1 the shock radii are plotted versus time for 
the upper and lower charges, and the two curves appear to be 
identical. Unfortunately, the relatively small charge 
separation in this experiment made it impossible to observe 
the primary shocks over a sufficient distance to calculate 
accurately the variation of shock strength with distance. 
The limited results which were obtained for shock strength 
versus distance are plotted in Figure 8.1, compared to 
Brode's (1957) calculations. 

3.3 Comparison of Mach shocks over different surfaces 

The refractive image analysis (Dewey et al, 1975) has 
shown what appears to be a significant difference between 
the strengths of the Mach shocks over the smooth ground 
and beneath the interaction plane between the two charges. 
The results of the particle trajectory analysis given in 
Figures 8.2 and 8.3 do not indicate the same difference. 
However, in the RIA case measurements were made as close 
as possible to the reflecting surfaces, 0.5m above the 
ground plane and 0.2m below the interaction plane, whereas 
in the PTA case the results represent an average of 
measurements made at puff positions which were not so close 
to those surfaces, at heights ranging between 1.0 and 4.5m. 


The results shown in Figures 8.2 and 8.3 therefore merely 


indicate that the difference in shock strength over the 


ground compared with that at the interaction plane may be 


dependent on the height above the ground at which the 
measurements are made - not an unexpected result. Determination 
of Mach shock strength from measurements made at various 
heights is also made difficult because an assumption must 
be made about the exact shape of the Mach shock front, in 
order to correctly assign shock radius values to smoke 
puffs in the PTA case. At or near a reflecting surface 
the problem of shape is not so important. Details of the 
problem in the PTA case and the manner in which the problem 
was dealt with for Shots 9 and 10 are described in Volume 1 
of this report. 
3.4 Resolution of time histories 

The time histories of density and pressure shown 
in Figures 24 and 25 do not always show a sharp rise at the 
time of shock front arrival. This slow rising is not a 
real effect but one inherent to the method of particle 
trajectory analysis, which does not permit a high 
resolution of density in space or in time. The finite spacing 
of the smoke puffs does not permit the average density of 
the air within a rectangular cell defined by four smoke 
puffs to be calculated accurately until the shock has completely 
traversed the cell. (The time of complete traversal may be 


“as much as 5 ms.) 


For the same reason the calculated time histories 
often anticipate the time of shock front arrival or, in 
other words, because a cell lies partly in front of the 
shock front during the time of traversal, the value of 
average density calculated at a point ahead of the shock 
may rise before the arrival of the shock. Also for 
the above reason, time histories calculated by the particle 
trajectory analysis method do not resolve any weaker 
shocks subsequent to the first, although these shocks 
may be detected occasionally in the calculated histories 
as a rounded bump in the normally exponentially decaying 
curve. Efforts are being made to see if this resolution 
can be improved. 

The lack of resolution close to the shock front does 
not occur in the case of particle velocity, which can be 
measured with reasonable accuracy as soon as the shock 
has traversed the relatively small space represented by an 
individual smoke puff. This effect of improved resolution is 
manifested also in the dynamic pressure histories since 
although this parameter depends on measured density, it 


also depends on particle velocity squared. The squared 


term is able to exert the greater influence. 
In the subsequent volume of this report hydrostatic 
overpressure and total head pressure histories determined 


from the particle trajectory analysis will be compared | 


directly with gauge measurements. The total head pressure 
will be calculated from the dynamic pressure using 
appropriate compressibility corrections. Pressure impulse 


curves will be included. 
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